Quantum Phases of Long Range 1-D Bose-Hubbard Model: Field Theoretic and 

DMRG Study at Different Densities 
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We use Abelian Bosonization and density matrix renormalization group method to study the effect 
of density on quantum phases of long range 1-D Bose-Hubbard model. We predict the existence of 
supersolid phase and also other quantum phases for this system. We have analyzed the role of long 
range interaction parameter on solitonic phase near half filling. We discuss the effect of dimerization 
in nearest neighbor hopping and interaction terms on the plateau phase at the half filling. 
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I. INTRODUCTION 

Different experimental and theoretical studies on su- 
perfluid and superconducting nano-scale systems reveal a 
rich quantum phase diagram (QPD) with many interest- 
ing quantum phases [H, H, H, 0, H, M, 0, [H] • One of the in- 
teresting quantum phases is the super solid (SS) phase in 
which the charge density and superconducting/superfluid 
phases characterized by diagonal and off-diagonal order 
coexist. The experimental findings and theoretical search 
for different quantum phases for cold atoms in optical 
lattice have revealed many interesting correlated phases 
of low dimensional bosonic systems [3, EH EH hj. In 
this regard, Bose-Hubbard model with extended range 
interactions have been studied in detail to discover the 
differe nt q uantum phases of cold atoms in optical lattices 
0, E3, El E2] • Here we study the quantum phases of a 
more general Bose-Hubbard model, namely, the Dimer- 
ized Bose-Hubbard model (DBH) with extended range 
interactions. The Hamiltonian of our model system is 
given by: 

H = -t 1 ^(l + (-l) i S t )(6tS i+1 +/ l .c) 

i 

- t 2^2(l>ih+2 + h.c) + — ^2hi(hi - 1) 

i i 
i i 

-fi^hi (1) 

i 

t\ and ti are the nearest-neighbor (NN) and next-nearest- 
neighbor (NNN) hopping terms respectively and V\ and 
V2 are NN and NNN interactions respectively. U is the 
on site repulsion energy and \i is the chemical potential. 
5 t and 5 V are the dimerization parameter for NN hop- 
ping and NN interaction respectively. Manipulation of 
interaction range and the prediction of different quan- 
tum phases in optical lattice loaded with cold atoms is 
more e asily constructed than other correlated systems 

[H M EE EE E3, El, El, M HH M M M Hi Dif- 
ferent combinations of laser beams with inhomogeneous 



intensity profile and their suitable manipulation can gen- 
erate long range interactions and anisotropic interactions 
extending to a desired range. So our theoretical model 
(DBH) is realizable because of the advances in the quan- 
tum state engineering of cold atoms in optical lattices. 
We believe that our theoretical prediction may help to 
understand and motivate experimentalist to design many 
interesting new systems. 



II. MODEL HAMILTONIAN AND 
CONTINUUM FIELD THEORETICAL STUDY 

Before presenting our numerical results, we briefly dis- 
cuss a field theory for the low energy and long wave length 
physics of DBH. We recast our basic Hamiltonian (Eq.l) 
in the spin language [26| to obtain; Hj 1 = —2 Ji 5^(1 + 

(-iys 1 )(s+sr +1 + h.c), Hj 2 = -2 J a 

h.c), H Ea , = E C0 E t ^S z l , H Eah = 4£ Zl £i(i + 

(-iys 2 )sf sf +1 , h Ec2 = 4Ez 2 J2i s i s i+2- 



H = H Jl + Hj 2 



Heco + Heci + He C2 



(2) 



The correspondence between the parameters of Eq. (1) 
and (2) is as follows: J\ ~ {n)t\ , J 2 ~ (n)t2 , 
Eco ~ (U < n > + //) , £7,i ~ V u E z2 ~ V 2 
[27| . One can transform the spin chain model to a spin- 
less fermion model through Jordan- Wigner transforma- 
tion with the relation between the spin operators and 
the spinless fermion creation and annihilation operators 
given by S* = iplip n -l/2 , 5,7 = tp n exp^X)?"^ n \ , 



St = tpl, exp[-iir Y% 



j] , [28| , where rij — ip^tpj is 



the fermion number at site j. We recast the spinless 
fermions operators in terms of field operators by the re- 
lation 

i>(x) = [e ik " x %l) R {x) + e* k " i, L {x)] (3) 

where iPr(x) and tPl{%) describe the second-quantized 
fields of right- (R) and left- (L) moving fermions respec- 
tively. We express the fermionic fields in terms of bosonic 
field by the relation 

U r 



1pr{x) 



-i (rcj>(x) - 6(x)) 



V2? 



(4) 



2 



r denotes the chirality of the R or L moving fermionic 
fields. The operators U r commute with the bosonic field 
as well as with U r of different species but anticommute 
with U r of the same species. (j> field corresponds to the 
quantum fluctuations (bosonic) of spin and is the dual 
field of (j>; 4>r = 9 — <j> and 4>l = + </>■ 

U sing the standard machinery of continuum field the- 
ory [28j |. we finally obtain the bosonized Hamiltonians 



Hi — Hn — 



S t Ji 
2na 



cos(2v Kcj){x)) dx 



(2na) 2 
AE Z1 5 V 



cos(4v K(j){x) ) dx 



(-If cos(4VK(t)(x) ) dx 



(5) 



where V{x) = £ r V r cos(^p), 



r, an integer and p[x) = 
Following Ref. 30 and 



31, the non oscillatory contribution of H2 arises from the 
commensurability condition nd — pa, d is the mean dis- 
tance between the particle, related to the density of the 
lattice. Under this condition, Hamiltonian for a particu- 
lar value of n is given by 



H2 = V n / dx cos(2p(j)(x)) 



(9) 



p = 1 is the most relevant commensurability and corre- 
sponds to one boson per site, p = 2 is the next relevant 
commensurability, with one boson every two sites. For 
these commensurabilities sine-Gordon coupling term be- 
comes relevant and system becomes gapped. 



''0 



Jo 27T 1 

2(^-2J 2 ) 



4(^1 - 2 J 2 ) 



dx : [DM*)? ■ + ■ [dxMx)f ■ 

Ji 



III. DMRG STUDY 

We now present numerical results obtained by using 
DMRG. We also compare them with the existing analyt- 
ical and numerical results. 



dx (d x (j) L (x))(d x (j} R {x)) 



(6) 



Hq is the gapless Tomonoga-Luttinger liquid part of the 
Hamiltonian with vo = sinfcp. The velocity, vo, of low 
energy excitations is one of the Luttinger liquid (LL) pa- 
rameters while K is the other. It reveals from Eq. 5 
that for weak dimerization, there is no contribution from 
the interaction part of Hi , given by the last term in Eq. 
5. The effective Hamiltonian obtained in this limit is 
the Hamiltonian for the saw tooth spin chain [2^ with 
dimerization. For strong dimerization, the Hamiltonian 
in Eq. [5] reduces to 



(2na) 2 



sia(4vK(j)(x) ) dx . 



(J) 



the second term in Eq. 5 and Eq. 7 yield a gap in the 
elementary excitations of system which led to plateaus 
in the fi vs p (boson density) in the system. In Density 
Matrix renormalization group (DMRG) study we will see 
evidences of plateau phases for different boson fillings and 
the effect of S t and 5y on these plateaus. We will also 
see occurrence of gapped phase for several commensurate 
fillings in our DMRG study, in the next section. Here we 
build up a general field theoretical study to explain the 
appearance of gap structure at different commensurate 
fillings: suppose we consider a periodic potential V(x) of 
periodicity of 'a' coupled to the density p{x) leading to 
an additional term in the Hamiltonian, 



H 2 = J dxV{x)p{x) 



(8) 



A. Numerical Details 

We use the Density Matrix renormalization group 
(DMRG) method to numerically study the QPD of the 
Hamiltonian in Eq. 1. We employ the infinite DMRG 
algorithm keeping 128 dominant density matrix eigen- 
vectors (DMEV) for determining p, while for the calcu- 
lation of correlation functions we use finite DMRG algo- 
rithm keeping the same cut-off in the number of DMEVs. 
Fock space of the site-boson is truncated to four states 
which allows 0, 1, 2 or 3 bosons per site. The length of 
the chain studied is 128 sites, except near phase bound- 
aries, where we have used 256 sites for calculating p 
and correlation functions. Accuracy of the method is 
checked by comparing the ground state (gs) energies, var- 
ious correlation functions and charge gap from DMRG 
studies with exact diagonalization studies of small sys- 
tems with upto 12 sites. We have also reproduced the 
results of earlier DMRG calculations satisfactorily [22| . 
The discarded density in the DMRG calculations is less 
than 10~ 14 in the charge density wave (CDW) phase 
at p = 0.5 as well as in the p = 1.0, Mott-insulating 
phase. However, the discarded density is slightly less 
than 10~ 10 in the superfluid (SF) phase. We have com- 
puted the charge gap (A) defined as A = limAr^oo A(AT); 
A(N) = E N (p) + E N (h) - 2E N (0), where N is number 
of sites on the chain and p, h correspond to an extra hole 
or extra particle at density p. E^(0) is the gs energy of 
zero particle or hole number at the same density. 
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FIG. 1: Energy gap (A) vs. density (p) for different values 
of Vi. The parameters for this figure are t\ — 0.1, t-2 = 0.0, 
U = 1.0, Vi = 0.7, 8 t = and S v = 



FIG. 2: Dependence of inverse correlation length l/£, S(7r) 
and |S , (2tt/3)| as a function of p for V2 = 0.0 and 0.45. Open 
circle corresponds to V2 = while square corresponds to 
1/2=0.45. Other parameter values are t\ = 0.1, U = 1.0, 
Vi = 0.7, all other parameters in Eq. [T]are set to zero. 



IV. RESULTS AND DISCUSSION 

Fig. Q] shows the variation of A with p for different 
values of Vi. At V2 = 0, we observe two peaks in the gap 
at the two densities, p = 0.5 and 1 as reported by Ba- 
trouni et.al [HI]. These peaks shift to w 0.35 and w 0.65 
on introducing nonzero V%. Position of peaks remains the 
same for the other nonzero value of Vi we have studied. 
We note from Fig. [1] that the gap occurs only near the 
two commensurate fillings of 1/3 and 2/3 (when Vi is in- 
cluded in the interaction) and disappears for fillings away 
from these values. This transition from gapped phase to 
gapless phase is the commensurate to incommensurate 
transition; the latter is due to the mismatch between the 
underlying periodic potential of the lattice and periodic- 
ity in the occupancy of the lattice. Fig. f2]shows the vari- 
ation of inverse correlation length of the density-density 
correlation functions l/£, as well as the structure factor 
S(q) computed for q = tt and 27r/3 as a function of boson 
density p, for two different Vi values, namely Vi = and 
Vi = 0.45. We first discuss the Vi = case. We note 
that for Vi = 0, the inverse correlation length shows a 
peak at p = 0.5 and 1.0 at which values we also note 
a gap in the system (Fig. [TJ. The underlying periodic- 
ity in the charge density at these p values corresponds 
to dimerization as seen from large S(tt) at these fillings. 
We also note that for Vi 7^ 0, the system has vanishing 
S(n) at all fillings. When Vi is switched on, the peaks in 



inverse correlation length shift to p = 1/3 and p = 2/3; 
at these values we also observe a nonzero gap (FigfT]) in 
the systems. The underlying charge order corresponds 
to a periodicity of three lattice sites for p = 2/3, as seen 
from the peak in \S{2n/3)\. We also note that \S(2w/3)\ 
is vanishingly small for all p in the case of Vi — 0. These 
results are also in broad agreement with field theoretic 
results. 

The charge charge correlation function in the gs for 
hole -doping and particle doping are shown Figs. [3Ja, and 
[SJa respectively. We note that for the case of hole doping, 
the 'defect' brakes up into two solitonic states each with 
charge half, for both values of Vi (0.1 and 0.3) for U = 1 
and V\ — 0.7. However, in case of the particle doping, 
we note that the two cases have quite different behavior. 
For Vi =0.1 we note that the charge-charge correlation 
function oscillates over the entire chain length. In case of 
Vi = 0.3, the oscillations are damped in the middle of the 
chain and become slightly more pronounced at the ends. 
This behavior is akin to what is seen in the hole-doping 
case. 

Physical picture for this behavior, can be arrived at 
from an analysis of the t = Hamiltonian. We find that 
at p = 0.5 the lowest energy configuration is the one in 
which alternate site are occupied by a single boson (Fig. 
2k top row) . On doping with a single hole we find that the 
energy reduces by 2Vi (Fig. fj^ middle row). However, if 
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j 

FIG. 3: Variation of density-density correlation with distance, 
(a) one less boson than half-filling, (b) one more boson than 
half-filling. Squares and circles represent V2 = and V2 = 0.3 
respectively. 
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FIG. 4: Schematic diagrams of different configurations of half 
filled bosonic chains with one extra particle and one extra hole 
together with their energies in the ti = £2 = limit. Panel 
(a) represents half-filling system with one extra hole whereas 
panel (b) is for one extra particle. 



the state with these consecutive holes is delocalized (Fig. 
12k bottom row) then there is a further stabilization by 
V-2- Thus the system prefers to break-up into two defects, 
with each defect corresponding to two consecutive hole 
sites. We can formally associate a charge half with each 
defect since two defects have been created by a single 
hole doping. 

The case of particle-doping is slightly different. When 
an empty site is doped (Fig. [Dd topline) then the energy 
increase is 2V\. Derealization of the particle, leads to a 
state with energy 2V\ — V2 ( Fig. [3J) middle line), which 
is stabilized by V% as in the hole-doped case. However, 
we can also dope a particle at a site which is already oc- 
cupied by a boson. The energy increase corresponds to 
U+ 2V2 in this case. Thus, we should observe a 3 consec- 
utive particle state yielding a state with two separated 
consecutive particle state ( Fig. |4Jd middle row) only for 
2Vi < U + 2V 2 - This is exactly what we find in Fig. 
HJx We can identify the ground state for 2Vi > U + 2V2 
as a Mott insulator state while that for 2V\ < U + 2V2 
corresponds to a solitonic state. We now turn our at- 
tention to the effect of dimerization on the phase dia- 
grams. We do not assume simultaneous dimerizations in 
both Vi and t\ as we wish to explore the role of each 
of these parameters independently. The effect of dimer- 
ization in t\ appears to smoothen the p vs. y, behavior. 
However, the dimerization in V\ seems to lead to higher 
jumps between plateaus, besides changing the value of 




P 



FIG. 5: Density (p) vs. chemical potential (/1) plot. Solid 
curve with circles are for nonzero dimerization in ti and 
dashed curves with squares are for nonzero dimerization in 
V\. Parameters space for this figure are t\ = 0.1, t-z = 0.0, 
U = 1.0, Vx = 0.7, V 2 = 0.45. 
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p at which the plateaus occur. It is easy to construct a 
real space picture for the observed p vs p behavior if we 
treat the transfer term as perturbation. In a model with 
next nearest neighbor interaction in the t\ = limit, it 
is possible to have a ground state with zero energy for 
all fillings, p < 1/3 for which the particle can be so dis- 
tributed that the inter particle interaction in Eq. [1] is 
zero. At, p = 1/3, when an extra particle is added, there 
is jump in the gs energy by (2V2) which is reflected as a 
step in the p vs. p plot. Further addition of particle will 
increase the gs energy by same amount until p < 1/2, 
keeping p constant between p — 1/3 and p = 1/2. This 
picture can be extended further for higher fillings. The 
effect of the transfer term is to reduce the sharpness of the 
jumps as well as introduce a slow variation in p between 
jumps. For p > 1/2, the chemical potential is nearly U 
for U < 2Vi and 2Vi for U > 2V\, when transfer term is 
switched on. Dimerization of the lattice does not signif- 
icantly affect the p vs. p behavior at least up to p = 1. 
Beyond p = 1, the analysis is not straightforward due to 
the large number of occupancy possibilities afforded by 
the the bosonic system. 

The physics of our system is similar to that of a sawtooth 
spin chain under a magnetic field mentioned in sec. II. 
Here p of the bosonic model is replaced by the magneti- 



zation and p is replaced by the magnetic field. We would 
like to give the physical explanation of the plateau state 
following the reference [32j ■ The energy levels of a mag- 
netic chain can be labeled by the M s value of the state. 
When an external magnetic field H is applied the state is 
stabilized by an energy —gHM s . Thus if the M s value of 
the gs in the absence of an external field is zero, when the 
field is turned on, the gs switches progressively to higher 
values of M s . If to begin with, the system had gaps 
between the lowest energy state in different M s sectors, 
then the M s value of the gs state shows jumps at discrete 
values of the magnetic field. This results in plateau in 
the M s vs H plot. 

In summary, We have carried out quantum phase analysis 
of dimerized Bose-Hubbard model, emphasing quantum 
field theoretic treatment as well as DMRG method to 
follow the quantum phases. A real space picture of the 
system in the zero hopping limit gives clear insights into 
the nature of the quantum phases, which are also pre- 
dicted by the quantum field in the strongly interacting 
limit. 
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